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We investigate the Levy glass, a mean-field spin glass model with power-law distributed couplings 
characterized by a divergent second moment. By combining extensively many small couplings with 
a spare random backbone of strong bonds the model is intermediate between the Sherrington- 
Kirkpatrick and the Viana-Bray model. A truncated version where couplings smaller than some 
threshold e are neglected can be studied within the cavity method developed for spin glasses on 
locally tree-like random graphs. By performing the limit e — > in a well-defined way we calculate 
£SJ ■ the thermodynamic functions within replica symmetry and determine the de Almeida- Thouless line 

in the presence of an external magnetic field. Contrary to previous findings we show that there 
is no replica-symmetric spin glass phase. Moreover we determine the leading corrections to the 
ground-state energy within one-step replica symmetry breaking. The effects due to the breaking 
of replica symmetry appear to be small in accordance with the intuitive picture that a few strong 
bonds per spin reduce the degree of frustration in the system. 
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INTRODUCTION 



Spin glasses have been investigated for about 40 years by now. The seminal analysis of Edwards and Anderson 
[l| revealed that disorder and frustration are the main ingredients necessary to bring about the peculiar static and 
dynamic properties of these systems. Subsequent analytical and numerical investigations have indeed shown that 
model systems with random interactions of the simplest possible type as, e.g., binary or Gaussian distributions may 
qualitatively reproduce various features of experimental spin glasses [H, Q . Trusting in universality and the ubiquitous 
efficiency of the central limit theorem no strong dependence of macroscopic properties on the details of the coupling 
O ■ distribution was expected. 

On the other hand many experimental realizations of spin glasses involve magnetic impurities placed at random in 
a non-magnetic metallic lattice. Mediated by the conduction electrons of the host material the impurities interact 
via the RKKY-interaction which oscillates in sign and falls off with distance r as 1/r 3 . In a homogeneous sample 
a given impurity hence interacts with order r 2 other impurities a distance r away. Correspondingly this impurity 
maintains order 1/ J 2 interactions of strength J. Probability distributions with such power law behavior are markedly 
\ different from simple ± J or Gaussian distributions. They describe a broad hierarchy of couplings and do not obey 
^vq ■ the central limit theorem Well-known representatives are Levy distributions characterized by a power law tail 
of the form l/J 1+a with the parameter a ranging between zero and two. Early investigations of spin glasses with 
RKKY-interaction p| Q already gave arguments for a broad distribution of magnetic exchange fields. More recently 
spin glass models with a wide hierarchy of coupling strengths have been used because they can lead to some models 
y—i , of finite dimensional spin glasses which may be well controlled (in the limit of a very strong hierarchy of couplings) 

> ; 0,1- 

• j-h ! The investigation of spin-glass models with power-law distributed couplings was initiated in 1993 by Cizeau and 
^ ■ Bouchaud H 10] who studied an infinite-range model using the replica symmetric (RS) cavity method. The model 
shows a transition from a paramagnetic high-temperature phase to a disordered glass phase at a freezing temperature 
T c depending on the parameter a. In their analysis of the low-temperature phase Cizeau and Bouchaud first provide 
arguments for a Gaussian distribution V(K) of local magnetic fields and then proceed to show that the model exhibits 
several unusual properties. The RS entropy becomes negative at sufficiently low temperature but returns to zero at 
zero temperature which gave rise to the speculation that RS may be restored at T = 0. Moreover, investigating the 
local stability of RS they found that the de Almeida-Thouless (AT) temperature Tat [H| is lower than the freezing 
temperature T c suggesting the existence of a finite temperature interval with a glass phase correctly described by RS. 

The model was re-investigated recently [L2l - [l4T | and it was found that the distribution of local fields is not a Gaussian 
and that the AT-temperature in zero field coincides with the freezing temperature excluding the possibility of an RS 
glass phase. In the present paper we give the detailed derivation of our results reported in [l4| and extend them in 
several directions. We provide a thorough analysis of the RS properties of the system, showing that the RS entropy 
does not vanish when T — > 0. We also characterize the correction which is to be expected from replica symmetry 
breaking (RSB) effects by determining the ground state energy of the system within one-step RSB. Most of our 
analysis is done in the framework of the cavity method, however, we make contact with the corresponding results 
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from the replica analysis. 

The Levy spin glass is intermediate between the two extreme prototypes of mean-field spin-glass models, the 
fully connected Sherrington- Kirkpatrick (SK) model [l5[ and the strongly diluted Viana-Bray (VB) model fla - flll ] 
respectively. Similar to the SK model in the Levy glass each spin interacts with all the other spins. The majority of 
the couplings is weak (©(TV -1 /")) as typical for fully-connected models. On the other hand due to the heavy tails 
of the coupling distribution on top of this background of weak couplings there is a backbone built from a few (0(1)) 
strong couplings per spin which remain 0(1) for N — > oo similar to the VB model. A somewhat related situation is 
given by composite systems |19| for which the properties of the weak and the strong bonds are defined separately. 

The decisive question is which properties of the Levy glass are exclusively determined by the strong bonds and 
which also feel the influence of the many weak ones. To elucidate this point we will often consider what we call the 
truncated model in which all couplings weaker than a certain threshold e are neglected. This technique has been crucial 
in the recent developments on the Levy spin glass problem, see [l3|, [l4| ■ We are then dealing with a spin-glass on an 
Erdos-Renyi random graph and use techniques developed for the cavity analysis of these systems [20L [2l| . Eventually, 
we have to investigate the crucial limit e — >• to recover the original Levy glass. 

The paper is organized as follows. In section II we introduce the model and the basic notation. Section III is 
devoted to the determination of the freezing temperature T c (a) and the RS distribution of local fields. In section IV 
we derive expressions for the thermodynamic functions like the free energy, the internal energy and the entropy and 
discuss the RS phase diagram of the Levy glass. Section V contains the determination of the AT-line of the Levy 
glass. In section VI we calculate the corrections to the ground state energy resulting from one-step RSB. Finally, 
section VII contains some conclusions. 



II. THE MODEL 



We consider a system of N Ising spins Si — ±1, i — 1, iV with Hamiltonian 

where the sum is over all pairs of spins, and /i ex t denotes an external magnetic field. The couplings Jij = Jji are 
independent, identically distributed random variables drawn from the distribution 

p a MJ) = ^jj^im-N-^), (2) 

where denotes the Hcavisidc function and the scaling of the couplings with N ensures that the free energy of the 
system is extensive. The most prominent feature of the coupling distribution ([5]) is its power-law tail for large values 
of | J | . In fact we will see that all macroscopic properties of the system depend only on the overall scale of couplings 
and the exponent a characterizing these tails. This implies in particular that all our result apply also to a spin glass 
with coupling distribution given by a symmetric Levy distribution defined by the characteristic function 

p}(q) = J dJPt{J)e^ J =exp(-^|C), (3) 

with 

~ 2sin(7ra/2)F(a + l) " ( ' 

In the large N limit the expectation values of any well-behaved function of J taken with respect to the distribution 
([2]) and ([3]) respectively coincide. 

In the present paper we will assume that a 2[ implying a finite average of \J\. With the Hamiltonian ([1]) being 
linear in the Jij we expect that for these values of a the thermodynamic potentials will be self-averaging. 



III. DISTRIBUTION OF LOCAL FIELDS 



A. Self-consistent equation for the local field distribution 



The central quantity in the replica symmetric cavity analysis of spin glasses is the distribution V(h) of local fields 
hi that parametrize the marginal thermal distribution of spin Si at site i 0, l2l| . Adding to a system of N spins 
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Si, i = 1, N another spin So with couplings Joi, i = 1, N one hnds for the local held at the new site the equation 

N 

h = h cxt + V u(hi, Joi) , (5) 



i=l 



where 



u(h, J) = ^artanh( tanh(^/i) tanh(^J)) . (6) 

The new field ho is a random quantity due to the randomness in the hi and Joi . The update equation ([5]) may hence 
be turned into a self-consistency condition by averaging over the distribution of hi and Joi and requiring that the 
statistical properties at site i = should be equivalent to those at all other sites. Accordingly 



V{h) = fYldh l V{h t )lY[dJo t P a ,N(Jo t )S(h~h cxt -J2u(h l ,Jo t )) 

i i i—1 

fdh'V(h')JdJP a . N (J) exp ( - isu(ti, J)) 



^■exp (is(h - h cxt )) 
Zn 



ds 
2^ 



JV-s-oo 



exp (is{h-h ext ) + N In 1 + ^ J dh'V{h') J ( cos ( s u (h' , J)) - 1)9 (| J\ - AT*) 

exp (is(h-h cst ) + | J dh'V{ti) J _^_(cob(« «(/»', J)) -1) 



(7) 



In the second line we used the statistical independence of the distributions at different sites. In the third one we 
inserted the explicit form ^ of P Qj at(J) and took advantage of the fact that it is an even function of J. Note that 
splitting off the leading 1 in the square brackets makes the J-integral well defined for N — > oo since cos(w( J, h)) — 1 is 
quadratic in J for small J and therefore suppresses the potential divergence of the integral for \J\ — > 0. This allows 
us to safely perform the limit N — > oo in the last line. 

Eq-0 is the central equation of the RS cavity approach. From its solution one gets V{K) which in turn determines 
all equilibrium properties of the system within replica symmetry. The same equation was also obtained in [l2| by 
using the replica method. 

For high temperatures and zero external field we must have (Si) = for all i and hence all local fields must vanish. 
Indeed V(h) = 5(h) is a solution of (|7|) for all temperatures. To test its stability one starts with a distribution V(h) 
with a small second moment and investigates whether it grows or shrinks under iteration. In this way one finds that 
the paramagnetic solution V(h) = 6(h) becomes unstable below the critical temperature 



T c (a) 



a / — tanh x 

X<*+1 



(8) 



This result was already obtained in Q and [Hj]. The dependence of T c on a is shown in Fig. [T] Note that T c diverges 
as a — t 2. On the other hand it is known that the Levy distribution tends to a Gaussian in this limit and we would 
hence expect to reproduce the results for the SK-model for a — > 2. This is indeed the case, however, the Levy 
distribution ([3]) approaches a Gaussian with divergent variance, J\, a =2, cf. (U). Rescaling T c with J\ a we indeed 
recover the classical results for the SK-model |15| . 

Below the freezing temperature Tn(a) no analytical solution of ([7]) is available. We note in particular that a Gaussian 
ansatz for V(h) as advocated in [9| does not reduce ([7]) to a self-consistent equation for the variance: Plugging in 
a Gaussian at the r.h.s. of ([7]) does not produce a Gaussian at the l.h.s. The only way to determine V(h) in the 
spin-glass phase T < T c (a) is hence by numerical methods. 

B. Numerical determination of the local field distribution 

We have implemented two ways to numerically solve ([Jj. The first consists in a straight iteration of the equation. 
From the n-th approximation V^ n \h) of the unknown distribution V(h) we determine 

Q(«)( s )=exp(| J dh'V^(h') J -^(cobM*',./))-!)) (9) 



4 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

a 



FIG. 1: Spin glass transition temperature T c (a) of an infinite-range spin-glass with couplings drawn from the distribution 
P a ,N defined in ((2J). 

by numerical integration. The next approximation, V^ n+ ^'[h), is then obtained via Fourier transform 

V {n+1) (h) = [ ^ e^-^xt) Q(»)( s ) . (10) 
J 2tt 

Starting initially from a uniform or normal distribution V^°\h) the procedure converges rather quickly. In order to 
save computation time we subdivide the J-integral in Q and approximate the small and large J part by analytical 
expressions. Moreover we use the Fast Fourier Transform with 256 or 512 nodes to perform the second step (fT0|) . 
Since in general results for several values of f3 are needed it is convenient to use the final result for one j3 as the 
initial distribution in the iteration for the next one. In this way smooth and accurate approximations for V(h) can 
be obtained in reasonable time. Moreover, in the limit j3 — > oo the integral in ([9]) may be simplified which makes the 
method very efficient for determining V{h) at zero temperature. 

Alternatively population dynamics as introduced for diluted spin glasses in [2(| may be used to solve ([7]). Compared 
to the direct iteration discussed above this method has two disadvantages: Firstly, it is statistical in nature and 
therefore one has to cope with intrinsic fluctuations. Secondly its application to the Levy glass needs the introduction 
of an additional cut-off parameter e for the coupling strength. On the other hand population dynamics has a big 
advantage which outweighs the above mentioned drawbacks: In a generalized form it may also be used for the 
investigation of the RSB phase (cf. section IVIj) . 

The most direct way to map the Levy glass onto a diluted spin glass amenable to population dynamics is by using 
the truncated model, i.e. by simply neglecting all bonds with modulus less than some threshold e. The number 
(K + 1) of remaining bonds per site is then a Poissonian random variable with mean . The distribution of the 
remaining bonds is given by 

^ )S (J) = ^^0(|J|- £ ). (11) 

Population dynamics may now be applied without further ado: Choosing K from its Poisson distribution and selecting 
at random (K + 1) values hk from an initial seed one replaces hx+i with 

K 

fencw = ^cxt + ^ u(h k , Jfe) (12) 
fc=l 

until the histogram of the hk no longer changes significantly. This procedure has to be performed for successively 
smaller values of e from which the asymptotic result for e — > may be extracted. 

Although this method works in principle, its convergence for e — > is slow. We found a significant speed-up of the 
algorithm by using the following modification. Instead of neglecting the weak bonds altogether we subsume them into 
a Gaussian random variable z with zero mean and a variance determined self-consistently. From the distribution 

p«M = N ' £ - a2lJ l +1 - N ~ 1/a ) ^ - \ J \) (is) 
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of weak bonds we find for this variance 



N- 



N \ 

u(hi,J 0i )\ =(N-e-°) I dhV(h) I dJP a , w (J)u 2 (h,J) 

= K+1 ) 

ihV(h)£^u 2 (h,J) =:^^L u \h,J)^. (14) 



Here as in the rest of the paper ( • ) h denotes the average with respect to the distribution of the local fields V(h) and 
the over-bar indicates the quenched average over the appropriate distribution of bonds. 
By combining the update (|14[) for z 2 with the noisy population dynamics algorithm 



A' 



J2u{h kl J k ) + z (15) 



fc=i 

we treat the dominant contributions from the strong couplings exactly, and include the influence of the weak couplings 
in an approximate way. This modification of (I12[) has at least two advantages. Most importantly its numerical 
implementation showed that the final extrapolation to e — >• is much smoother. Moreover, in the opposite limit, 
e — > oo, all couplings are included in the Gaussian variable z and we may check the algorithm by comparison with 
the results obtained in [9|. 

The integrals entering (|14[) are still computer demanding. To save some computer time it is useful to tabulate the 
function 

a dJ 2 



H ^ Jo J^ u (h ' J) 

for an appropriate interval of values of the local field h before the update procedure. The determination of the variance 
([T4"f is then reduced to a simple integral at each update. 

The left part of Fig. [5] shows the distribution of local fields V(K) in zero external field as obtained by direct iteration 
of eq. (JT)). We find practically indistinguishable results by running 10 3 iteration of population dynamics with e < 0.3 
using 10 4 members in the population. Markedly different, however, are the Gaussian distributions proposed in Q 
which are also shown in Fig. [5J The reason for these differences lies in the fact that that the sum in ([5]) is dominated 
by a few large contributions. Consequently, although the different terms in the sum are statistically independent and 
all have finite second moments, the central limit theorem may not be invoked since the Lindeberg criterion is not 
fulfilled. 

To further check our numerical results we have determined the second moment ofV(h) as a function of temperature 
close to T c (a), cf. the right part of Fig. O The variance tends to zero when the temperature approaches T c from 
below as it should. Moreover, the slope coincides with the one following from the analytical expansion of ([7]) for small 
T — T c which gives 

(h 2 ) h = ~T e {a) (T - T c (a)) + 0((T T c {a) f) . (16) 

Our solutions for V(h) are similar to those given in Fig. 2 of [Hj]. However, our results for the Gaussian approxi- 
mation are significantly different from those shown there. 



IV. REPLICA SYMMETRIC THERMODYNAMICS 



In this section we give expressions for the thermodynamic potentials - free energy, internal energy and entropy - 
which, on the replica symmetric level, are all functionals of the local field distribution V{h) determined above. We 
first derive a differential equation for the free energy per spin f(j3, /i ox t) of the model defined by eqs. (JJ) and ^ within 
the cavity method. Next we consider the truncated model and derive an expression for the free energy and the internal 
energy using the cavity method for diluted spin glasses. We then show that, in the limit e — > 0, this free energy fulfills 
the differential equation derived before. We also show how the same expressions for the thermodynamic functions can 
be derived from the replica method. Finally we discuss the most salient features of the RS thermodynamics of the 
Levy spin glass. 



h 



T 



FIG. 2: Left: V(h) for a = 1.1, T = 0.9 T c , 0.6 T c , 0.1 T c (center from top to bottomV Full line: results from the iterative 
solution of the self-consistent equation (|7). Dotted line: Gaussian approximation from 9]. Right: Second moment of V(h) for 
a = 1.1, 1.5, 1.8 (from left to right) as a function of temperature. Symbols: As determined from the iteration results for V(h). 
Lines: From the analytical expansion of (O around T c , cf. (|16l) . 



A. Differential equation for the free energy per spin 



Assuming that the free energy per spin /(/3, /i e xt) is self-averaging in the thermodynamic limit it can be related to 
the shift in free energy due to the addition of one spin via 



log Z N (J3, ftext) 



Pf{P,K*) = ~ lim b JV ^' =- lim (logZ N+1 ((3,h ext ) - logZ N (P,h ext ) , (17) 

where the disorder average is taken with respect to the coupling distributions P a> N+i and P a ,N respectively. Due 
to the explicit dependence of the coupling distribution on the number of spins the comparison between systems of 
different size needs some care Q. In the present case the slight change of the coupling distribution ^ when going 
from N to (N + 1) spins is absorbed in the rescaled parameters 



N j -ext y N 

Splitting off the terms depending on the new spin So we find 

The clustering property of the pure state in the absence of the new spin, i.e. the statistical independence of all spins 
interacting with the newcomer implies 

™=n™> (i9) 

which when used in (|18|) yields 

N „ N 



log ^ff^y = ^logcosh(/3J 0l ) + ^log(l-tanh 2 (/3^)tanh 2 (/3J 0l )) + log (2cosh(/3/i )). 
Averaging (|20|) over the disorder then gives 



log Z z + N \f h ^f = N J dJ ^(J)logcosh(/3J) + <log(2cosh(/3/ l ))) h 

N f 

+ Y J dJ P Q ^ N {J) (\og(l - tanh 2 (l3h)tanh 2 ((3J))) h . (20) 
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The limit N — > oo may now be taken on both sides of this equation. Using (fT7|) and neglecting 0(1/ N) contributions 
we find for the l.h.s. 

l0 S ry ,n 7 T~ = ~pf + oY ~ = ~Pf ( e + ^ext™) , 

Z N {p,h ex t) a dp a oh cxt a 

where we have used the thermodynamic relations for the internal energy and magnetization per spin respectively 

d(pf) R d(pf) 

e = \ a P m = — -■ 

dp dh cxt 

Combining this result with the limit of the r.h.s. gives rise to the differential equation 



pf+^(e + h cxt m) = -(log(2cosh(^))) fc -/ 



adJ 

2T7F+1 



logcosh(^J) + -/log (l - tanh 2 (/3/i) tanh 2 (/3J)) 

2 \ / h 



.(21) 



For zero external field this equation was already derived in |l0| where, exploiting the assumption of a Gaussian 
distribution of local fields, also an explicit solution was constructed. However, in view of the fact that the correct 
distribution of local fields is non-Gaussian and is not available analytically, a straight integration of (f2"Tj) to find 
/(/?, h ex t) is difficult. Instead we will use two different ways to derive an expression for the free energy per spin and 
verify that it indeed fulfills (I2T1) . 



B. Free energy of the truncated model 

When neglecting all bonds with a strength less than a threshold e the Levy spin glass is converted into a spin glass 
on a locally tree- like random graph Q N e - a with N sites, mean connectivity £~ a , and a coupling distribution given 
by (fTTj) . We may therefore use methods from the cavity analysis of the Bethe spin glass [2(| with only minor changes 
due to the fluctuating connectivity in our model. The free energy per spin of the truncated model is given in terms 
of two different free energy shifts according to 



1 . t- = : 1 



Pfe(h ext ,p) = -^ lim (logZ N+2 -logZ N ) = £ 1(K + 1)AfP -2KAF™) . (22) 

Here AF^ corresponds to the free energy shift due to the addition of a single spin and AF^ to that due to the 
addition of two spins connected by a bond. An intuitive explanation for this relation can be obtained by considering 
two operations acting on the graph Q N e - a of the truncated model. Removing 2K vertices from this graph leads to a 
cavity graph where some spins lack neighbors, here as before (K + 1) is a poissonian with mean e~ a . Adding (K + 1) 
new pairs of spins ao, tq connected by a bond Jo to the system, and connecting them to the free spins produced by the 
first operation leads to a Gn+2,£-" graph where the connectivity remains unchanged whereas the number of vertices 
is increased by 2. The resulting expression for the free energy per spin reads : 

J dJP a , s {J) logcosh/3J+y J dJP a , s (J) (log(l -tanh 2 (/3J)tanh 2 (/3/i))) h 



(log 2 cosh ph) h - — — I d.JP a , s {J) (log (l - tanh 2 (/3J) t&nh 2 (ph) tanh 2 (/3/i')))^ h , 



Taking into account the £-dependence of the distributions of the connectivity and the couplings strength the limit 
e — > may be performed and we obtain for the free energy per spin of the original model 

M = -J I ^^[logcosh(/3J) + (log(l-tanh 2 (/3J)tanh 2 (^))) /i ] -(log(2cosh(/3/ l ))), i 

(log (1 - tanh 2 (/3J) tanh 2 (/^) tanh 2 ^')))^, ■ (23) 

We now turn to the determination of the internal energy e £ per spin. There are two contributions: one due to the 
interactions and one due to the external field, e E = ej, mk — h cxt m. To obtain the link contribution we follow the steps 
in (20| and add a coupling to the system. For the Levy case it is convenient to rewrite the expression for the 
energy of this link obtained in [2(| as 



Jij (SiSj) =-Qjp log ( cosh(/3' J y )(l + tanh(/3%-) tanh(/3/iJ tanh(/3^))) 
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here hj and hj denote the local fields in the absence of the new link at the site i and j respectively. A link connects 
two spins, each of which interacts on average with e~ a = K + 1 neighbors. After the average over the quenched 
disorder the system is homogeneous. The link contribution to the internal energy is hence related to the average of 
Eij by 



Jink 



K + 1 



2 2 

1 r adJ d 

'2j\j ]>e 2\J\^W 



- [ dJP a , s (J) (^log(cosh(/3'J)(l + tanh(/3'J)tanh(/3/i)tanh(/3/i'))) ) 

J \VP ^ 'I h.h' 



13=13' 



logcosh(/3'J) + - (log (1 - tanh 2 (/3'J) tanh 2 {/3h) tanh 2 (/3h'))) h h , 



13=13' 



where we skipped all contributions which vanish due to the symmetry of the P a _ s distribution. Using in addition 
m = (tanh(/3/i))t we find for the internal energy of our original model : 



e = limel ink 



h cxt (tanh(/3/i)) h 



dp 

a 
'2P 



7?" 





f adJ 


I3=fl' 2 J 


2|J|«+! 



logcosh(J) + - (log (1 - tanh 2 (J) tanh 2 (f3h) tanh 2 (0h'))) h 



h cxt (tanh(/3/i)) ?l 



a dJ 
2|J|«+ 1 



logcosh(^J) + - (log (l - tanh 2 (/3J) tanh 2 (/3/i) tanh 2 {/3ti))) h h , 



h ext (tanh(0h)) h . (24) 



The results obtained for /, e and m fulfill the differential equation (|2"Tj) derived in IIV Al The replica symmetric 
thermodynamics of the truncated model is hence in the limit e — > equivalent to that of the Levy spin glass. Note 
that the above reasoning relies on the fact that the limits TV — > oo and e — > commute. However, since all expressions 
depend smoothly on the cutoff parameter for e — > we believe that this is indeed the case. In order to further 
substantiate our results (|2"5|) and (IM1) we re-derive them in the next subsection using the replica approach. 



C. Free energy from the replica approach 



Within the replica approach the free energy per spin is related to the replicated partition function Z n via 

13 f = - lim Jim — log Z n . 



(25) 



Due to the slow decay of P a ,N(J) the quenched average of Z n does not exist for n ^ 0. As in [12| we therefore use 
imaginary temperatures ft = —ik , k G R at intermediate steps of the calculation, i.e. before the limit n — > is taken. 
For integer values of n, the quenched average of the partition function reads 



Z n {-ik) = f Y[ dJij P a , N (Jij) exp(-ikJSi ■ Sj) 



(26) 



J2 ex p U^ 108 '! 1 



{S°} 



idJ 



N I 2 J 



cos (kJSi-Sj) - 1 



0(N~<) 



where 7 < 1 for all values of the parameter a considered here. Eq. (|26|) can now be transformed into a 2™-dimensional 
integral over order parameters [22| 



1 N 



(27) 



where a — {& a } denotes an n-component Ising vector and S is the Kronecker-5. The partition function acquires the 
form 



Zn = J ]Jdc(a)d(^2c(a)~?j exp (- N(-ik) j trial ({c (a)})) , 



(28) 
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and is evaluated by the saddle-point method for N — > oo. The free energy per spin is then determined by 

/= ]im-f tria i({co(a)}) (29) 

where Co minimizes the trial free energy 

Hfc)/ 4rwi ({ C (a)})=^c(a)logc(a)-i^c(a) C (a') J [cos ( k J ff ■ ff') - 1] . (30) 

The corresponding saddle-point equation reads 

= l + A(n)+log Co (?) + ^co(5 ! ') J [cos ( k J ff ■ ff') - 1] (31) 

where the Lagrange multiplier A(n) accounts for the constraint c(ff) — 1 resulting from (|27|) . Within the replica 
symmetric assumption Co(ff) depends on the sum Y^l=i a a only, and is related to the distribution V(h) of local fields 
via [12 

c« s (ff) = / AhV{h) ^- lkh ^\ a : ] . (32) 
J (2cosh(-ikh)) 1 ; 

The saddle-point equation ([31]) for Co(ff) can then be transformed into a self-consistent equation for V(h) which 
coincides with Q. Using the RS ansatz ([32l the sums in (f30|) can be performed leading to expressions for which the 
n — > limit can be taken (for details see [12J ). The value of the Lagrange multiplier A(n) can be inferred from the 
saddle-point equation at a a = 0. In the limit n — > we find 

lim 1 + ^ (n) = (log2cosh(/3fo)) ft + J logcosh^ + i^l-tanh^^tanh 2 ^)))^ (33) 

where the continuation to real temperatures has already been performed. Using ([29|) . (|3T1) . and ([32)) we find back 
expression (1231) for the free energy per spin. To obtain the internal energy we use 



1 1 f) I \ 

e ^ ^iV^^^^o^^ E^-p UE« ■ (34) 

{5 ? } V o=l / 

Again the limit n — > and the continuation to the real temperatures may be performed and the result for the internal 
energy obtained in section ITVBI gets reproduced. 



D. The RS thermodynamic functions 



In the previous subsections the thermodynamic functions of the Levy spin glass within the assumption of replica 
symmetry were determined using different approaches. In the left part of Fig[3]we have plotted the free energy per 
spin for zero external field, h ex t = 0, as function of temperature for three different values of a. In the right part of 
this figure the internal energy per spin in zero external field and at T = is shown as function of a. For small T 
the numerically obtained values for the free energy smoothly approach those for e(T = 0). The RS free energy has 
negative slope for high temperatures, reaches a maximum, and has a negative slope near T = 0. Correspondingly 
the entropy becomes negative at sufficiently low temperature, a well-known signature for the breakdown of replica 
symmetry in spin glasses. We have plotted the entropy per spin in the inset of Fig. [3] Instead of using the derivative 
of /(/?) with respect to T it is numerically much more accurate to determine the entropy from the thermodynamic 
relation s — j3{e — /). For temperatures close to T c (a) the curves for the entropy are similar to those obtained within 
the Gaussian ansatz for the distribution of local fields. At lower values of T, however, there are significant deviations 
from the results obtained in |9J. In particular we neither find a minimum of s(T) at low temperatures nor do our data 
extrapolate to s(T = 0) = 0. 

The temperature T s=n at which the entropy becomes negative decreases with decreasing a, cf. Fig. [3] This is 
consistent with intuition since a larger fraction of strong bonds should reduce the degree of frustration. We will find 
a similar behavior when studying the influence of RSB in section I VII Moreover we find for all a that T s= q is much 
smaller than T c . In I9j it was speculated that the strong bonds in a Levy glass may stabilize a RS glass phase for some 
finite temperature interval below T c . To investigate this question we study the stability of RS in the next section. 



10 



-2 




FIG. 3: Thermodynamic functions of the Levy glass as calculated from the numerical determination of the order parameters 
using population dynamics. Symbols give error bars of the statistical error intrinsic to population dynamics, lines are guide for 
the eye. Left: Replica symmetric free energy per spin for a Levy glass with a = 1.8, 1.5, 1.1 (green, blue, red). To lighten the 
comparison the data have been normalized to the ground state energy f(T — 0). The inset shows the corresponding results for 
the entropy per spin. As characteristic for spin glasses the replica symmetric entropy becomes negative at low temperature. 
Right: Replica symmetric ground-state energy per spin, e = f(T — 0), for a Levy spin glass as function of a. The red symbols 
show the ground state energy for a = 1.8, 1.5, 1.1 within one-step RSB. As can be seen the corrections are rather small (cf. 
also Fig.©. 



V. STABILITY OF THE RS SOLUTION 



The self-consistency of the RS cavity approach can be tested by investigating the correlations between spins. More 
precisely, the divergence of the spin-glass susceptibility 

^G = ^E(^>-^>^->) 2 (35) 

(i,3) 

signals the breakdown of replica symmetry Q . In order to determine the stability boundary for the Levy spin glass we 
start again with the truncated model for which all weak bonds smaller than e are neglected. We then use techniques 
from the theory of diluted spin glasses [2l[ to determine the region of validity of RS and finally perform the e — » 
limit. As we will see this limit may be accomplished analytically which makes the extrapolation back to the original 
model safe. 

The underlying graph of the truncated model is locally a tree which allows to write (|35l) as 

oo 

XSG = Y, £ ~ ar Cr, (36) 

r=l 

where 

C r = (<5 5 r )-(So>(5r)) 2 (37) 

denotes the square of the connected correlation function of two spins at a distance r and e~ ar gives the average 
number of sites at distance r from i — 0. For large r we expect C r ~ exp (— r/£) with some correlation length £ 
and therefore the divergence of the sum (f5rjf depends on whether £ -a e -1 '^ ^ 1. The stability analysis for the fully 
connected Levy glass is thus mapped on the asymptotic behavior of the spin glass correlation in a one-dimensional 
disordered Ising chain, cf. Fig. @] It is defined by the Hamiltonian [2lj | 

r — 1 r 

E({Si}) = -J2 J&Si+i - hiSi ' ( 38 ) 

1=0 8=0 

where the couplings Ji are distributed according to the distribution P a>s of strong bonds (fTTj) and the external fields 
hi are drawn independently from the distribution V(h) determined in section fill Bl 
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FIG. 4: Left: Part of the graph corresponding to the truncated model. The couplings Ji on the red branch are treated exactly, 
whereas the influence of all other couplings on the spins of this chain is approximated by the local fields hi sampled from V(h). 
Right: The corresponding one-dimensional model defined by (|38|l . 



The correlations in the one-dimensional system defined by (|38| have been studied already in [23] using transfer 
matrix techniques for the replicated systems. We re-investigate the problem by using a cavity approach which leads 
to the same result but sheds light on the importance of rare fluctuations. Let us start again with an update equation 
of the form for the local fields of model 



9i+i = hi+i + u(gi, Ji) . (39) 



According to linear response theory the correlation function Ci is related to the change of magnetization at site i due 
to a perturbation 5h of the local field at site 0. Denoting by ( • ) S the canonical average induced by the Hamiltonian 

E s ({Si}) = E({Si}) - 5hS (40) 



we hence have 

i 2 

.9 g 



'd(Si) 5 



dSh 



/3 2 [1 - tanh 2 (/3 5l )] 2 ( ) =: 0' cosh' 4 ^) D l (41 ! 



5/1=0/ 

Using (|39l) and the chain rule we can derive a recursion relation for Di : 

A + i=(^|^) 2 A- (42) 

The quantities gi and Di are hence correlated random variables due to their dependence on the quenched disorder 
represented by the fields hi and the couplings Ji. From the update rules (|3"9"|) and (|42l) we find for their (site dependent) 
probability distributions 

P i+1 (g i+1 ,D i+1 ) = J dJiP a , s (Ji) J dh i+1 V(h i+1 ) J dgt J dD l P l (g l ,D l ) 

5(g i+x - [hi+x+uiguJi)]) s(D i+x - A)- (43) 

This equation describes how a perturbation propagates through the chain, and thus contains information on how the 
correlation C r = cosh _4 (/3g J ,) D r decays with increasing r in a given sample. It could easily be studied with population 
dynamics in order to get the typical decay rate. However, in order to characterize the spin-glass susceptibility P6|) . 
we need the behavior of the average of C ri not the typical one, and this average correlation is dominated by rare 
instances of the couplings. 

In order to obtain the behavior of the average correlation we first derive a recursion relation for the auxiliary 
quantity 



Ii(gi) := / dAPi(ft,A) A (44) 
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from which the averaged correlation function can be obtained by integration, C r = J dg cosh 4 ((3g) I r (g)- The average 
decay of correlations is hence determined by the i-dependence of Ii(g). We now find from (|43j) 



I i+ i(g i+ i) = j dJiP a , s {Ji) J dh i+1 V(h i+ i) J dgi J dDj Pi(gi, A) 



/ dujgj, Jj) 

V % 



A S {g l+1 - [h i+1 + u(g i: Jj)]) 



/ du(g u J) 



J d 9i J dj p a , s (j) J dhv(h) ih: 



s {di+i ~[h + u(g l , J)}) Ii(gi) 



(45) 



from which we infer that the asymptotic behavior of h{g) and hence also that of the averaged correlation function 
C r for r — > oo is characterized by the largest eigenvalue v of the transfer matrix K(x,y) defined in the last step of 
(|45]) . We therefore conclude that the stability of RS is determined by the convergence of the geometric series 



XSG 



E 



r = l 



(46) 



The most convenient way to determine the stability parameter A is to calculate the largest right eigenvalue of e a K, 
i.e. to solve the equation 



\cj>{x) 



dy K T (x,y) cj>(y) 



tdJ 



2\J\ C 



dhV{h) 



du(x, J) 
dx 



<fi(h + u(x, J)) 



(47) 



where K T denotes the transposed operator, and the e — > limit has been taken in the last expression (notice that it 
can be taken safely since the function (du(x, J)/dx) 2 behaves as J 2 for small \J\). For A < 1 RS is stable, otherwise 
it is unstable. The correlations in the one-dimensional system defined by (|38[) have been studied already in 23] using 
transfer-matrix techniques for the replicated systems. The eigenvalue problem for the replicon mode considered there 
coincides with (|47p. 

The determination of the largest eigenvalue in f|47[) has to be done numerically and can be accomplished by straight 
iteration. Starting with an arbitrary positive function <f)^ we use 



a dJ 
2\J\ a +- 



dhV{h) 



( du[x, J) 
\ dx 



^(m-l) (h-\- U (x, J)), 



(48) 



and impose the normalization J 0( m ) (x) dx = 1 after each step. After many iterations </>( m ) converges to the eigenvector 
of the operator (ITT)) with the largest eigenvalue which in turn is given by A = linim^oo Z m . The complete numerical 
procedure is hence as follows. For given f3 = 1/T and h ex t we calculate the distribution of local fields V(h) from Q 
as described in section ITlI B I We use V(h) to determine the transfer matrix K(x,y) according to ([4"5"j) and extract its 
largest eigenvalue X(T,h ext ) from the iteration (|4*5|) . The AT-line is then implicitly given by A(/i cxt , Tat) = 1- 

In the left part of Fig. [5] we have plotted the stability parameter A as function of T — 1/(3 for different values of h cxt . 
One finds that A(T, ft ex t) is, for all values of h cxt , a monotonically decreasing function of temperature and crosses the 
line A = 1 at exactly one point. The collection of these points defines the AT-line TAT(h cx t) shown in the right part 
of Fig. [5] for three values of a. 

Of particular interest is the value of Tat in zero external field. In the paramagnetic region characterized by 
V(h) = 5(h) the stability parameter may be determined analytically. Using u(Q, J) — we find in this case from (|47|) 
for x = 



A0(O) 



a dJ 



dhS(h) 



( du(x, J) 
\ dx 



4>{h + u(0, J)) 



x=0 



a dJ 



T t'&n\i{l3J)<i){Q) 



T 



0(0), (49) 



Assuming 0(0) ^ we hence find the instability of the paramagnetic solution at T c (a). Moreover, our numerical 
results indicate that A(T, 0) > 1 for all T < T c (a). 

Because the numerics is somewhat subtle when V(h) is near to a 5-function we corroborate this result by a pertur- 
bative study of the eigenvalue problem to leading order in the reduced temperature r := 1 — T/T c (a). To this end 
we expand (|47p in x up to order x 4 . Using the symmetry of V and cj> the truncated eigenvalue equation acquires the 
form X(f> — K</> with the vector <p = (0(0), <fi"{0), <fi""(0)) T . The matrix elements Kij depend on the moments of V 
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FIG. 5: Determination of the de Almeida- Thouless line. Left: Stability parameter A as function of temperature for a = 1.5 and 
foext = 0, 0.5, 1 and 2 (from right to left). From the intersection of the curves with the stability boundary A = 1 the AT- line is 
determined. Right: Phase diagram of a Levy spin glass with a — 1.8, 1.5, 1.1 (from right to left). Above the AT-lines shown 
RS is stable, below it is unstable. 



which to the required order in the reduced temperature read 



/ h 2\ a rj.2, \ , a rp2 f \ a(ll*2,a + 28t 4iQ ,) - 18(t 2) , 



■ T 



V /h t 2 , a -U,a V 2^ cV ' 



t2,a — tA, 



(50) 



where we defined 



For the largest eigenvalue of K we find 



t m , a ■= tanh m ( J) . 



A = 1 



+ 0{t a ) 



(51) 



(52) 



The coefficient of the quadratic term is always positive, since t m , a decreases with m as implied by (1511) and tanh(x) < 1. 

We therefore find A(T, h ex t = 0) > 1 for all T < T c (a) in the Levy spin glass. Correspondingly there is no stable 
replica symmetric glass phase as proposed in [l(| ■ The AT- line for the Levy spin glass as shown in the right part 
of Fig. [5] looks indeed qualitatively similar to other spin-glass models. 

We finally comment on the assumption cf>(0) ^ made after eq. (|49|) . For <p(0) = we differentiate (|47|) m times 
where d m (j> is the first derivative with d m cf)(0) ^ 0. Evaluating this equation at x — in the paramagnetic region i.e. 
for V(h) = 8(h) we infer that the corresponding eigenvalue A TO is given by 



Am, — 



a dJ 



tanh 



2+m 



t(2+m),, 



Since t m . a decreases with m, it results that eigenfunctions with (f)(0) = give rise to eigenvalues which are smaller 
than those for eigenfunctions with <p(0) ^ and are hence not relevant for the stability problem at hand. 



VI. ONE-STEP RSB 

The RS solution fails at low temperatures as it ignores the possibility of several pure states [Iq| . We therefore 
consider a solution which takes into account the existence of many pure states and corresponds to the one-step 
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RSB solution in the replica formalism. The local fields h] corresponding to the different states 7 are assumed to 
be independent random variables sampled from a site dependent distribution Pj. After averaging over the disorder, 
the natural order parameter is the probability distribution Q of the local field distributions Pi, and one can derive 
a self-consistent equation for it. In this section we work at zero temperature, where RSB effects should be most 
pronounced, and we show how to compute the ground-state energy density of the model within the one-step RSB. 
For simplicity we shall keep to the case of zero external field. The method which we use to derive expressions for the 
ground-state energy consists in using the truncated model, to which the general one-step RSB approach to diluted 
models developed in [24j can be applied. The limit e — > is performed analytically at the end of the calculations. 
We shall derive the complete one-step RSB equations and solve them numerically using an approximation called the 
factorized approximation [25| . 

The assumption of one-step RSB is the following: for a fixed realization of the disorder the local fields h] on a given 
site i are random variables due to the existence of many local ground states. We denote by Pi the corresponding site 
dependent probability distribution. 

When iterating, i.e. merging K spins at a new site, the update rule for the probability distribution at this new site 
reads 



K K K 

PnewiKew) = C / dhi P t (hi) 5 (h ne w ~ ^u (Jj, hi) J exp (/iy^max(| Jj\ - \hi\,0) + fJ,\h(,\), 



(53) 



where (i is the one-step RSB parameter, the constant C ensures the normalization of P neu ,, and uqQi, J) is the zero 
temperature limit of ([5]) u (h, J) — lim^_ i . 00 u(J, h) = min(|/i|, |J|) . (Notice that we use a slightly different notation 
from the one in [24j : we consider the update rules for the probability distributions P%(hi) = CjPj(/ij) exp(— fi\hi\) , 
instead of Pi(hi) used in [24[; this is to ensure a safe e — >■ limit). 

The ground-state energy of the Levy spin glass is obtained within the one-step RSB approximation by maximizing 



$( / u) = lim $ £ ([i) = lim 



AE° ite (p) - -(K + l)AE b e ond {^i) 



(54) 



with respect to fx, where the energy shifts A£'| Jte (/x) and /S.E bond (fi) are given below. 
The energy shift corresponding to a site addition reads 

AEf e ({Pi,Ji}\fi) = log / Y[dh i Pi(hi)expUj2 max (\ J i\-\hi\,0)+^ 

^ J i=l \ i=l 



K+l 



u {Ji,hi 



(55) 



i r K+i 

-logC-log J dH P {H )e-^ - ^log J dhiPi(hi) exp(nmax(\ Ji\ - \hi\,0)) , 



where Pq denotes the distribution of the fields Hq = ^fj^ uo(Ji, hi) acting on the new spin connected to (K + 1) 
old ones. 

The energy shift due to a bond addition amounts to 



AE b E ond (P 1 ,P 2 ,J\fi) = --log / dhdtiP 1 (h)P 2 (ti)ex-p(LLm&x(h<7 + tia' + JcTcj')-n\h\~Li\ti 



= -|J|--log 



dhdti Pi(h) P 2 (ti) [l + 6»(-J/ l / l ')(e- 2 ^ min (l J l^ i l'l ,l 'l) -1) 



(56) 



After the average over disorder, the added site has the same properties as the old ones, and a self-consistency equation 
for the order parameter Q can be derived: in ([53)) P new must have the same distribution Q as the {Pi}. Assuming 
that Pq is also from the same distribution Q (which is the case in the e — > limit) the averaged energy shifts read 



AEf e {ii) = i (logJdhP(h) 



e -»\h\ _ £ -c 



P a:S (J)log(l + J dhP(h) [epv-MW-Wfi) _ i]^ (57) 



and 



AE>> '(/,) = / (1J 1],.aJ)( |J| + ^log / dhdti P(h)P'(h') \ +6(-Jhh')(e 2 " 



(58) 



p.p' 



respectively, where ( • )p denotes an average in which P is drawn from the probability distribution Q. 
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We have now all ingredients entering $ £ (/i) of the truncated model and the limit e —> can be performed to obtain 
the corresponding expression for the Levy spin glass. The symmetry of the averaged field distributions is crucial in 
this case as 

l^J dhdh'P(h)P'(h')9(-hh')j =~ (59) 
leads to a convergent integral in (|58|) at small values of J. One obtains 



oK\J\-W) 



$( M ) = i^log J dhP(h)e-»W- J™ ^-log^l + J' AhP{h)[e 

l^{^2 + h l0g l ^ P{h)P ' {h,) V + 0(-hh')(e- 2 ^^ h ^ 1)] ^> 




p,P' 



Notice that the self-consistent equation for the mean P with respect to Q, denoted by P(h) = J dP Q[P] P(h), 
coincides, in the limits e — > and fi — > 0, with the self-consistent equation ([7]) for the RS order parameter V for zero 
temperature and zero external field. In particular one recovers in this limit: 

lim $0) = lim fusifi, h ext = 0) (60) 

ft— >0 13— >oo 

A. Factorized approximation 

It is numerically rather heavy to sample P from the distribution Q[P]. The factorized approximation makes the 
task much easier by confining the space of all distributions to one distribution P, i.e. it assumes that Q is a functional 
delta-function. Within this ansatz a population dynamics algorithm can be easily applied to determine P for all 
values of the parameter fi. One then has to find the maximum of the function 

1 r roc i t 

S/GO = -log / dhP(h)e-^ + / _ RAJ) (61) 



where i?„ p denotes 



R,A J ) = -^ lo g( 1 + fdhP(h)[e^ J -W-l])+^ 



+ —log [ dhdti P(h)P(ti) \l + 6(-hti)(e- 2f * mia( - J ^^ h '^ - 1 
2fj, J L 

An expansion of i?„ p at small and large values of J allows to perform the J integral in (paT]) analytically in the 
corresponding regions. For small values of J one obtains R p(J) ~ \^J 2 , whereas in the large J limit one has 

R„,p( J ) ~ -4-- lo § l dhP(h)e-^ + ^-\og(l + \ [ dhdh' P(h)P(h')[e- 2 ^ min( - w ^ h ' l) - 1] 

A numerical test of the result for large J shows that the asymptotic behavior is reached for J > 50. Using this 
decomposition, the computation of ( E>/(/^) becomes numerically easy. We have obtained the distribution P by a 
population dynamics algorithm after performing 10 4 iterations for a population of 10 5 fields. The result is shown 
in Fig. [6l which plots the relative increase compared to the RS solution as a function of the one-step RSB 

parameter fi. The maximal value of $f(n) gives the estimate for the ground state energy of the Levy spin glass within 
the one-step RSB factorized approximation. We see that the relative corrections to this ground state energy due to 
RSB remain rather small, of the order of one to two per cent, which is comparable to the typical corrections found 
for instance in the SK model [26j . It is also smaller when a is close to 1, which agrees with the intuition according to 
which a smaller value of a leads to a stronger hierarchy of couplings and therefore to a lesser degree of frustration. 
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FIG. 6: Trial ground energy $/(a0 of a Levy glass as function of the one-step RSB parameter /i for a — 1.1 (red), 1.5 (blue), 
and 1.8 (green) as obtained within the factorized approximation. Symbols show population dynamics results with error bars, 
lines are guides for the eye. The curves are normalized to the modulus of the RS ground state energy = 0)|. The ground 

state energy in one-step RSB is given by the maxima of the curves which is about one per cent higher than the RS result. 



VII. CONCLUSION 

We have presented a detailed study of the properties of the Levy spin glass at the replica symmetric level at all 
temperatures and magnetic fields, and a one-step RSB study at zero temperature and zero external field. One main 
ingredient of this study has been the introduction of a truncated model where the couplings with values smaller than 
a cutoff e are neglected. The truncated model naturally enters the category of dilute spin glasses for which various 
techniques have been developed in recent years, allowing for a detailed analysis. The e — > limit requires some care, 
and complicates notably the analysis with respect to the studies of 'usual' spin glasses, but we have shown that it can 
be controlled. 

The physical picture which has been obtained shows a spin glass behavior which is generally much closer to 
the standard behavior found in the SK or in diluted models than what had been claimed before. Within the RS 
approximation, the entropy decreases with temperature and becomes negative at low temperatures, but does not turn 
back to when T — > 0. The AT instability line can be computed, and the whole spin glass phase turns out to be 
unstable with respect to RSB effects. On the other hand, the quantitative effects of RSB on the ground state energy 
are relatively small and become smaller with decreasing the Levy exponent a. 
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